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Abstract 


A  model  was  developed  for  determining  the  space  and  time  behavior  of  the  magnetic  field 
external  to  a  system  of  moving  pulsed  finite  conductivity  conductors.  The  moving  conductor  is 
in  relative  close  proximity  to  a  stationary  conductor.  The  external  magnetic  field  content  has  a 
space  and  time  periodicity  that  is  related  to  the  parameters  of  the  system.  A  numerical  example 
predicts  the  magnitude  and  frequency  of  the  magnetic  field  as  a  function  of  distance  from  the 
stationary  conductor.  Potential  applications  of  this  study  are  discussed. 
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1.  Introduction 


Over  the  years,  there  has  been  continuing  interest  in  the  development  of  pulsed-power 
energy  sources  and  accelerators  based  on  the  processes  of  transient  magnetic  field 
exclusion  and  diffusion.  An  essential  feature  of  these  devices  involves  a  stationary 
source  current  and  the  generation  of  induced  currents  (i.e.,  eddy  currents)  in  a  moving 
finite  conductivity  conductor.  The  source  and  eddy  currents  produce  magnetic  fields 
external  to  their  geometry  that  have  a  characteristic  magnitude  and  frequency. 
Applications  utilizing  this  process  can  be  found  in  electromagnetic  braking  of  large 
electromechanical  systems  [1,  2,  3]  and  the  generation  of  high  electrical  power 
pulses  [4, 5].  These  applications  require  100' s  of  kA  with  millisecond  pulse-widths. 
Additionally,  these  currents  can  produce  external  magnetic  fields  of  a  few  Tesla. 

As  these  electrical  devices  mature  into  components  for  military  systems,  we  must 
consider  the  electromagnetic  compatibility  (EMC)  issues  that  arise  from  the 
“environmental”  electric  and  magnetic  fields  they  generate.  By  “environmental”  fields, 
we  mean  those  fields  that  are  created  in  spatial  locations  that  are  not  especially  close  to 
the  regions  that  they  originate  from  and  are  concurrently  located  in  close  proximity  to 
electronic  equipment  that  may  be  affected  by  them.  For  example,  we  have  found  that 
environmental  electromagnetic  fields  generated  by  capacitor-driven  pulsed-power  sources 
could  be  important  EMC  considerations  for  electromagnetic  railguns  [6,  7].  Fortunately, 
standard  shielding  techniques,  using  conducting  materials,  are  adequate  for  managing  all 
but  the  largest  magnetic  fields  (i.e.,  <  IT)  [8].  For  large  confined  fields,  active  shielding 
techniques  can  be  more  efficiently  implemented  [8]. 

The  purpose  of  this  initial  assessment  is  to  characterize  the  salient  features  of  the 
environmental  magnetic  field  generated  by  a  moving  conductor  system  that  is  produced 
by  a  spatially  periodic  time-varying  source  current.  Figure  1  shows  a  close-up  of  the 
relationship  between  the  spatially  periodic  source  currents  and  the  moving  conductor. 
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Figure  1.  Source  Current  Elements  and  Moving  Conductor. 


Figure  2  shows  the  geometry  used  for  determining  the  environmental  magnetic  field. 
The  object  of  this  study  is  to  determine  the  time  behavior  of  the  magnetic  field  at 
point due  to  the  source  currents.  In  Figure  2,  the  system  is  assumed  to  be 

periodic  and  infinite  in  the  ±  jc -directions.  We  have  cast  the  problem  in  cartesian 
coordinates  for  mathematical  simplicity,  but,  with  the  proper  mathematical 
transformation,  it  can  be  made  to  correspond  to  cylindrical  geometry. 


In  this  report,  no  attempt  is  made  to  relate  the  test  results  to  a  specific  application.  A 
significant  amount  of  detailed  information  would  be  required  to  accomplish  this,  and  such 
an  effort  is  beyond  the  scope  of  the  present  study. 

The  results  of  our  model  quantify  the  damped  harmonic  electromagnetic  fields 
generated  by  a  moving  conductor  in  terms  of  the  periodic  spatial  separation  between  the 
current  sources,  the  conductivity  and  velocity  of  the  moving  conductor,  and  the  location 
of  the  observation  point  relative  to  the  moving  conductor.  The  principal  contributions  to 
the  far-field  magnetic  signature,  however,  are  shown  to  be  due  to  the  source  current 
distribution  and  its  time-dependent  image  created  during  a  transient  pulse. 
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Figure  2.  Geometric  Considerations  for  Determining  the  Environmental  Magnetic 
Field  Produced  by  a  Moving  Conductor  System. 


The  mathematical  details  leading  to  the  determination  of  the  current  density  in  the 
moving  conductor  and  the  subsequent  prediction  of  the  environmental  magnetic  field  are 
tedious  to  follow,  hi  order  to  facilitate  the  process  of  understanding  the  basic  physical 
ideas,  analysis,  and  interpretation  of  the  results,  a  brief  overview  of  the  basic  physical 
processes  is  presented  in  section  2.  The  theoretical  prediction  of  the  current  density  and 
associated  magnetic  fields  in  the  moving  conductor  is  rendered  in  section  3.  Using  the 
results  of  section  3,  the  environmental  magnetic  field  is  computed  in  section  4. 
Numerical  results  are  rendered  in  section  5.  Summary  and  conclusions  are  presented  in 
section  6. 
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2.  Theoretical  Considerations 


Much  insight  into  the  structure  of  the  electromagnetic  fields  described  in  Figures  1 
and  2  can  be  gleaned  by  first  considering  the  solution  for  a  single-source  current  element. 
For  convenience,  we  select  the  center  source  element  of  Figures  1  and  2.  Elucidation  of 
the  essential  physical  phenomena  is  further  facilitated  by  considering  the  two-dimensional 
case  in  which  source  current  elements  are  infinite  in  the  ^-direction  and  where  the 
thickness,  A ,  is  suitably  large.  The  implications  of  this  latter  assumption  are  defined  in 
the  course  of  the  discussion. 

For  illustration,  we  begin  with  the  case  in  which  the  conductor’s  velocity  is  zero.  We 
want  to  determine  the  response  of  this  nonmoving  conductor  to  a  sudden  increase  in  the 
source  current,  /(?)  :t>0,  which  flows  through  each  source  element.  (The  origin  of  I(t) 
is  not  essential  for  understanding  the  issues  related  to  the  electromagnetic  fields. 
Nevertheless,  it  is  worthwhile  to  mention  that,  in  the  case  of  an  energy-producing  device, 
l(t)  is  determined  from  external  considerations  that  involve;  generation  of  an 
electromotive  force  in  the  stationary  conductor,  an  electrical  load,  and  the  current  induced 
in  the  moving  conductor  by  I(t)  itself.  For  example,  in  the  case  of  a  compensated  pulsed 
alternator  [5],  the  effect  of  the  moving  conductor  is  to  produce  a  time-dependent  mutual 
inductance.  Detailed  discussions  on  the  types  of  compensation  and  the  methods  for  the 
physical  realization  of  compensation  in  pulsed  machinery  can  be  found  in  the 
literature  [5].) 

The  problem  of  the  stationary  conductor  further  divides  into  two  categories:  (1)  the 
case  where  the  conductivity  of  the  moving  conductor  is  infinite  and  (2)  the  case  where  the 
conductivity  is  finite.  A  close-up  of  the  situation  for  a  single  conductor  with  a  =  <»  is 
shown  in  Figure  3. 
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Figure  3.  Single  Source  Current  Element  and  Its  Im^e  for  an  Infinite  Conductive 
Surface. 


The  basic  physical  law  that  determines  the  current  flow  in  the  case  where  a  =  oois 
concerned  not  with  the  magnetic  field,  but  with  the  electric  field.  The  surface  impedance 
for  a  medium  that  has  infinite  conductivity  is  identically  zero.  This  requires  that  the  total 
component  of  the  electric  field  on  the  surface,  which  is  parallel  to  the  surface,  must 

also  be  identically  zero.  A  reflected  electromagnetic  wave  must  therefore  have  its 
component  of  the  surface  electric  field  that  is  parallel  to  the  surface  exactly  cancel  the 
corresponding  component  of  the  incident  electric  field  [9].  The  only  component  of  the 
electric  field  that  can  exist  on  the  surface  of  a  perfect  conductor  is  the  component  that  is 
normal  to  the  surface,  . 


The  foregoing  discussion  concerning  applies  at  all  frequencies.  It  can  be  shown 
using  the  theory  of  a  horizontal  dipole  over  a  conducting  half  space,  taken  in  the  limit  of 
<7  =  oo ,  that  the  condition  E^=0  leads  to  image  theory  [9].  When  applied  to  Figure  3, 
image  theory  shows  that  the  magnetic  field  at  the  surface  of  an  infinite  conducting  body 
can  be  computed  from  two  sources  [9].  The  first  source  is  the  hue  physical  current 
source  (flowing  out  of  the  paper),  and  the  other  is  the  image  of  the  source  current  that 
flows  into  the  paper. 


5 


Image  theory  permits  the  field  at  all  points  in  space  to  be  computed  by  summing  the 
fields  due  to  the  actual  physical  real  source  and  its  image.  The  real  source  and  its  image 
each  produce  a  magnetic  field  determined  from  the  Biot-Savart  law  [9].  The  three 
dimensional  form  of  the  Biot-Savart  is  given  by 


Ux,y,z)xr 

-5 - wcdyazt 


(1) 


where  X  ts  the  current  density;  x,  y,z  are  the  coordinates  of  the  volumetric  distribution 
of  X ;  and  X,  y,  z  are  the  coordinates  at  which  the  field  is  to  be  computed; 

r  =  T(x-x)  +  j(y-y)  +  k(z-z)  (2) 

and 

r  f  =  ((;r  -  x)^  -I-  (y  -  yf  +  (z  -  zff^ ,  (3) 

where  T ,  j ,  and  k  are  unit  vectors  in  the  x-,  y-,  and  z-directions,  respectively.  For 
currents  that  flow  only  in  the  z-direction,  we  have 

7c  =  Jc(x,y,z)k  (4) 

and 

kxr  =  j(x-x)-T(y-y).  (5) 

When  applied  to  the  two-dimensional  configuration  of  Figure  3,  the  contributions 
from  the  real  source,  ,  and  image,  Bj ,  are  given  by 
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(6) 


+/j  /  2  +0O 


B. 


-/,/2-«.((x-x)  +(y  +  w)  +(z  — z)  / 


and 


+/j  /  2  +O0 


]  It — 

_,,/2-i((jJ-x)  +(y-w)^+(z-z)  J 


■dzdx. 


(7) 


where 


7.(0  = 


/(O 

h 


(8) 


is  the  current  per  unit  length  of  the  source  current  element.  When  equations  (6)  and  (7) 
are  evaluated  on  the  surface  of  the  conductor,  }?  =  0 ,  we  find  that  the  ^’-component  of  the 
magnetic  field  vanishes  and  the  ^-component,  ,  consists  of  equal  contributions  from 
the  real  source  and  its  image.  We  get 


/  2  +0O 


B, 


(x)  =  -^J„{t)w  I  j- - — - ^ 

2^  -fnio^x-x)  +(wf +(z-z)J 


■dzdx. 


(9) 


Integration  over  z  gives 


+1,12 


^  I  / 


_l,^{x-xf  + 


w 


^B(-x). 


(10) 


Using  the  principle  of  superposition,  the  total  surface  magnetic  field  corresponding  to  the 
alternating  array  of  Figure  1  is 
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(11) 


n=-Ko  r 

n=:-oo  ^ 

By  making  the  substitution  x  =  jr  +  mL ,  where  m  is  any  integer,  in  equation  (11),  we 
easily  see  that  is  a  periodic  function  of  L.  We  derive 


^  n'L, 

BHx  +  mL)=  ZB„(x-—X-ir' 


(12) 


where  the  new  dummy  variable  of  summation  is  now  n'  —  2m- n.  Since  n  and  n'  are 
summation  indices,  equation  (1 1)  is  the  same  as  equation  (12). 

Combining  the  synmietry  and  periodic  properties  of  ,  we  can  write 


~  'ZtITC  X 

BL  =  Ssi(Ocos- 

n=l 


L  ’ 


(13) 


where 


B"  =  —  \ 

«  T  J 


2  f  2n7Cx 


cos- 


B^dx. 


(14) 


The  modal  components,  ,  are  determined  by  straightforward  integration  of 
equation  (14)  using  equations  (10)  and  (11).  The  details  of  the  calculation  are  not 
important.  It  is  important  to  note,  however,  that  all  are  proportional  to  the  product  of 
-  .  Accordingly,  it  is  convenient  to  use  the  form 

(15) 
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where  the  Q„  are  shown  to  be  dimensionless  expansion  coefficients  given  by 


The  magnetic  field  intensity  at  the  surface  due  to  all  the  conductors  is 


^  Innx 

K=^=  COS— — 


(17) 


B=1 


The  general  expression  for  the  surface  current  density,  7, ,  is  given  by  the  vector  equation 

nxH^=J„  (18) 

where  is  magnetic  field  intensity  on  the  surface  and  n  is  the  normal  to  the  surface  of 
the  conductor.  As  shown  in  Figure  3,  n  points  in  the  negative  y-direction  and 

71=00  OwTT  T 

H,  =  iHl  =  -r/„(t)Sn.  cos—— .  (19) 

n=l  ^ 


Inserting  equation  (19)  into  (18)  gives 

7  Tt  /  'IflKX 

Js  =  COS——  . 

n=\  ^ 


(20) 


In  contrast  to  the  discontinuous  elements  of  the  array  shown  in  Figure  1,  the  current 
density  on  the  surface  is  always  continuous  because  the  surface  is  continuous.  For 
example,  in  the  case  of  a  single  element,  which  is  discontinuous,  the  surface  magnetic 
field  has  been  shown  in  equation  (10)  to  be  continuous.  Using  equation  (10)  in 
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equation  (18)  then  provides  a  continuous  surface  current.  Discontinuities  in  surface 
current  lead  to  singular  point  accretions  of  charge  on  the  surface.  This  does  not  appear  to 
be  physically  plausible  in  this  problem. 

Equation  (20)  shows  that  J ^  is  a  periodic  function  of  x.  By  comparison,  the 
discontinuous  current  density,  J  ,  of  the  periodic  distribution  of  source  elements  shown 
in  Figures  1  and  2  can  be  written  as  =  kJj^ .  The  magnitude  is 

4(^»0  =  SA.(0cos-— ,  (21) 

n=0  ^ 


where 


L  =  2(/,  +  4), 


(22) 


and  the  \  ’s  are  determined  from  the  orthogonality  properties  of  the  Fouier  series.  Using 
the  discontinuous  distributions 


7(a:,0  =  /„(0 


and 
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(23) 


we  derive 


A.= 


2/.(0 

nn 


(l-cosn;r)sin 


nizli 

~ir' 


(24) 


The  image  distribution  for  the  array  of  source  elements  is  the  negative  of  equation  (21). 


In  contrast  to  the  image  current,  which  is  a  mathematical  abstraction,  the  surface 
current  of  equation  (20)  is  a  real  current.  However,  from  the  point  of  view  of  calculating 
the  conductor’s  contribution  to  the  magnetic  field  for  y  <  0  from  the  Biot-Savart  law,  we 
have  a  choice  of  using  either  the  real  surface  current  located  at  y  =  0  or  the  image  current 
located  in  the  plane  y  =  w. 

In  summary,  we  have  learned  the  following  from  the  infinite  conductivity  case  that 
has  bearing  on  the  general  case  involving  a  moving  conductor  having  finite  conductivity. 

(1)  A  true  physical  surface  current  is  produced  on  the  conductor  that  is  generated  by 
the  source  current  array. 

(2)  The  physical  surface  current  is  computed  from  the  source  current  and  its  image. 

(3)  The  surface  current  is  continuous. 

(4)  The  time  behavior  of  the  surface  current  and  its  image  is  proportional  to  the  time 
behavior  of  the  source  current.  For  example,  if  the  source  current  were  to  remain 
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constant  in  time,  so  would  the  surface  current  and  the  image  current.  The  reason 
for  this  is  that  infinite  conductivity  precludes  energy  dissipation  and,  hence,  a 
current  once  generated  will  persist  indefinitely  (in  the  absence  of  radiation  losses). 

We  now  consider  a  nonmoving  conductor  with  finite  conductivity.  This  extension  is 
one  that  involves  an  enormous  increase  in  complexity  as  compared  with  the  infinite 
conductivity  case.  As  previously  mentioned,  the  basic  physical  problem  to  be  solved  here 
is  that  of  a  horizontal  dipole  over  a  conducting  ground  plane.  This  problem  was  first 
attacked  by  the  German  theoretical  physicist,  Arnold  Sommerfeld,  at  the  beginning  of  the 
20*  century,  and,  for  the  first  half  of  this  century,  played  a  crucial  role  in  the  field  of  radio 
science.  Analytic  study  of  the  famous  “Sommerfeld  Integrals”  connected  with  the 
solution  continues  to  this  very  day.  The  seminal  text  by  Banos  [10]  provides  a 
comprehensive  discussion  of  this  fundamental  and  difficult  problem. 

The  solution  for  the  horizontal  dipole  over  a  conducting  ground  furnishes  a  complete 
space-and-time  description  of  the  electric  and  magnetic  fields  both  inside  the  conductor 
and  in  free  space.  For  example.  Hill  and  Wait  [11]  solve  explicitly  for  the  fields  in  the 
case  where  the  line  source  lies  directly  on  the  ground.  Our  problem  is  simply  an 
extension  of  the  basic  problem  to  include  a  finite  number  of  distributed  elements. 
Although  this  type  of  rigorous  formulation  can  be  attempted,  it  is  not  worthwhile  to  do  so 
unless  we  can  demonstrate  that  simpler  models  can’t  do  the  job. 

The  literature  on  horizontal  dipoles  over  ground  planes  is  extensive  and  has  been 
addressed  by  the  leading  electromagnetic  theoreticians  in  this  century.  A  major  outcome 
of  these  studies  has  been  the  emergence  of  an  approximate  image  theory  that  is  applicable 
for  conducting  media  [12].  Figure  4  shows  how  this  theory  is  applied.  If  w  is  the 
location  of  the  image  for  <7  =  <» ,  the  total  displacement,  ,  for  a  «>  is  given  by  [12] 

w^  =  w->rd,  (25) 
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Figure  4.  Approximate  Image  Theory  for  Conducting  Ground  Plane. 

where  d  is  the  additional  displacement  due  to  the  conductivity.  In  frequency  space,  d  is 
actually  a  complex  quantity,  a  consequence  of  fitting  the  rigorous  theory  with  a 
mathematical  approximation  [12].  It  is  given  by 


2 


1/2  > 


(26) 


which  is  valid  for  the  condition,  a  »(oe  (applicable  in  our  case). 


If  we  associate  the  magnitude  of  d  with  the  physical  displacement  and  let  (O^ln/  tp, 
where  is  the  duration  of  a  pulse,  we  obtain 


w= 


(27) 


Equation  (27)  is  (within  a  numerical  constant  of  order  unity)  the  root-mean-square 
displacement  determined  from  diffusion  theory.  If  is  small,  the  displacement  of  the 

image  will  likewise  be  small.  Equation  (27)  thus  provides  the  physically  intuitive  result 
that,  for  short  times,  the  location  of  the  image  is  that  obtained  from  the  <7  =  oo  case.  In 
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order  to  apply  equation  (27)  with  any  degree  of  confidence,  the  thickness,  A ,  must  be 
greater  than  \d\. 

Li  order  to  use  image  theory,  we  need  to  know  not  only  its  displacement,  but  also  the 
magnitude  of  the  image  current.  Whereas  in  the  <t  =  oo  case,  the  magnitude  of  the  image 
current  is  equal  to  that  of  the  source  current,  this  will  not  rigorously  be  true  in  the  finite 
conductivity  case.  The  magnitude  of  the  image  current  will  be  less  than  that  for  the 
infinite  conducting  case  because  there  are  energy  losses  in  the  system  associated  with 
finite  conductivity.  The  rigorous  solution  of  this  problem  requires  the  formulation  in 
terms  of  a  frequency-dependent  horizontal  dipole  over  a  conducting  half  space  [10]. 

When  the  velocity  is  not  zero,  the  location  of  the  image  is  the  same  as  that  for  the 
stationary  conductor  for  a  =  <» .  The  reason  for  this  is  the  same  as  before:  the  only  way 
we  can  assure  zero  electric  field  at  the  surface  is  for  the  image  to  be  opposite  to  the  real 
source.  Since  the  source  is  stationary,  so  is  the  image. 

An  explanation  of  current  behavior  when  the  conductor  is  moving  and  a  0°  is  not 
easily  explained  without  solving  the  entire  problem.  Conceptually,  this  problem  can  be 
attacked  as  a  moving  horizontal  dipole  over  a  conducting  half  space.  To  our  knowledge, 
this  approach  has  not  been  attempted.  By  all  reckoning,  it  would  appear  to  be  an 
extremely  difficult  problem.  Nevertheless,  this  is  the  problem  that  is  addressed  in  this 
study. 

Our  solution  is  based  on  solving  Maxwell’s  equations, 

V^H=aE  +  oVxB  ,  (28a) 

VxE—,  (28b) 

dt 
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(28c) 


and 


V*5=0  (28d) 

in  a  conductor  moving  with  velocity,  V .  In  our  model,  it  is  assumed  that  the  conductor’s 
thickness,  A ,  is  large  enough  so  that  the  electromagnetic  fields  never  diffuse  out  of  the 
moving  conductor.  Mathematically,  the  large  enough  assumption  means  that  we  can  set 
the  magnetic  field  components  equal  to  zero  at  y  =  A. 

The  foregoing  assumption  about  A  being  large  has  only  an  indirect  bearing  on  the 
calculation  of  the  environmental  magnetic  field  from  the  currents  in  the  system.  What  is 
really  important  is  the  depth  of  penetration  of  the  current  into  the  moving  conductor  as 
compared  with  the  vertical  distance  of  the  field  point  at  which  the  magnetic  field  is  being 
calculated  from  the  conductor  surface.  As  we  show  in  diis  study,  the  depth  of  magnetic 
field  penetration  is  quite  small  compared  with  meaningful  distances. 

An  accurate  prediction  of  the  fields  is  possible,  provided  that  we  can  precisely  specify 
the  magnetic  field  at  the  bottom  surface  for  all  time,  f  >  0 .  This  is  exactly  where  the 
difficulty  arises;  the  magnetic  field  at  y  =  0  cannot  precisely  be  defined  for  all  time 
without  simultaneously  solving  Maxwell’s  equations  in  the  space  below  the  moving 
conductor. 

It  would  appear  that  an  approximate  solution  to  this  problem  can  be  obtained  using 
perturbation  theory.  The  concept  is  as  follows:  first,  we  solve  for  the  electromagnetic 
fields  in  the  conductor  by  assuming  that  the  magnetic  field  at  y  =  0  is  given  by  image 
theory  for  a  =  «> .  This  approximation  enables  us  to  determine  the  first-order  correction 
for  the  displacement  of  the  image  due  to  diffusion  and  the  decrease  in  the  magnitude  of 
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the  image  current  due  to  energy  losses  associated  with  finite  conductivity.  Corrections  to 
the  surface  magnetic  field  are  then  made. 


3.  Electromagnetic  Fields  and  Currents  in  a  Moving  Conductor 


Maxwell’s  equations  in  a  moving  conductor  are  given  by  equation  (28).  In  the 
geometry  of  Figures  1  and  2,  the  velocity  is  given  by  the  relationship:  V  =  i  v .  By  taking 
the  curl  of  equation  (28a)  and  making  the  substitutions  of  equations  (28b)  and  (28c)  for 
VxE  and  H  in  equation  (28a),  we  derive  the  following  identical  equations  for  and 


and 


1 

lia 


dB^ 

- —+  V - - 

dt  dx 


(29a) 


1 

fio 


dB  dB 

V - L 

dt  dx' 


(29b) 


Equations  (29a  -  b)  are,  however,  not  entirely  independent.  The  divergence  condition 
of  equation  (28d), 


V*B=0=%+^=0,  (30) 

dx  dy 

creates  an  interdependency  on  the  solutions  of  equations  (29a  -  b).  They  are  each  second- 
order  in  X  and  y,  and  first-order  in  time.  If  it  were  not  for  the  divergence  relationship, 
they  would  be  completely  independent.  For  the  boundary  conditions  of  our  problem,  we 
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readily  show  that  the  solution  of  equation  (29a)  combined  with  equation  (30)  uniquely 
determines  By. 

Let  us  first  consider  the  solution  for  .  We  must  specify  the  boundary  conditions  at 
y  =  0  and  at  y  =  A .  The  boundary  condition  at  y  =  A  is  greatly  simplified  for  pulse 
lengths,  tp,  that  are  short  (i.e.,  A  >1  I).  When  these  conditions  apply,  we  are  at  liberty 

to  set  both  components  of  the  magnetic  field  equal  to  zero  at  the  top  boundary: 

B^{x,y  =  L,t)  =  By{x,y  =  h,t)  =  0.  (31) 

This  is  the  boundary  condition  used  in  our  study. 

The  foregoing  boundary  will  be  relevant  only  for  those  cases  where  the  diffusion 
time,  from  the  initiation  of  the  pulse  at  y  =  0  to  y  =  A  is  larger  than  tp.  The 

diffusion  time  is  approximated  by  the  formula 

(32) 

Diffusion  time  can  vary  greatly,  depending  on  the  application.  For  example,  using  an 
aluminum  conductor  (<7=2.2x10’  mhosim,  /i„=4;rxl0“’  H/m),  and  a  nominal  thickness 
A=  1.5  cm  gives 


tj)=6.2ms.  (33) 

When  tp  is  significantly  shorter  than  t^,  the  top  boundary  of  Figure  1  is  never 

reached  during  its  duration  and,  for  all  practical  purpose,  it  could  be  located  at  y =<» .  On 
the  other  hand,  for  a  system  in  which  there  is  a  continuous  stream  of  pulses,  the  condition 
at  the  top  boundary  may  need  to  be  more  critically  examined. 
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The  solution  of  the  problem  requires  that  we  specify  the  time  behavior  for  at  y=0 
for  r  >  0 .  This  quantity  is  denoted  in  this  section  as 

B°  =  B^{x,y  =  Q,t).  (34) 

B°  is  determined  from  the  considerations  rendered  in  section  2.  For  example,  we  can  set 

2,n7tx 

COS——,  (35) 

(1=1  ^ 

which  is  the  field  due  to  the  periodic  source  of  Figure  1  over  a  perfectly  conducting  plane. 
Use  of  equation  (35)  constitutes  the  first  approximation  in  a  perturbation  scheme. 

Our  analysis,  however,  is  carried  out  in  general  terms  by  keeping  the  spatial 
periodicity  and  allowing  for  an  arbitrary  time  behavior  in  7^(r) .  This  is  accomplished 
using  a  Delta  (impulse)  function  current  density  source.  The  boundary  condition  at  y  =  0 
then  becomes 

n=^  r 

cos—— .  (36) 

n=l  ^ 

If  B^  and  By  are  now  interpreted  as  the  solutions  to  equations  (29a)  and  (29b)  for  the 
impulse  function,  the  actual  fields,  B^(x,y,t)  and  By{x,y,t)  will  be  given  by 

t 

B^{x,y,t)  =  J  B^{x,y,x)J„(t  -  x)dT ,  (37a) 

0 

and 
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(37b) 


By{x,y,t)  =  j  By{x,y,T)J„{t  -  x)dx . 

0 

In  section  4,  we  apply  the  foregoing  formulas  to  the  case  where  Jgit)  is  a  step  function 
that  has  the  properties 


J^  =  Aj  :t^>t>0  (38a) 

J,=0  (38b) 


where  Aj  is  a  constant. 


Denoting  the  Laplace  transform  variable  by  “s”  gives 


1  ^Bx^ 

B(P\_  ^  By^  \ 


=  sBx  +  V 


(39) 


where  Bx  is  defined  as  the  Laplace  transform  of  Bx.  The  foregoing  equation  is  next 
reduced  to  one  involving  only  the  y-dimension  by  using  the  periodicity  of  the  x-direction. 
We  have 


BJ,x  +  L,y,s)  =B^{xy,s) . 


(40) 


Using  equation  (40),  we  next  write  Bx  as  the  real  part  of  the  harmonic  series 


JTnX 


(41) 


where 
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(42) 


Substituting  equation  (41)  into  equation  (39)  and  equating  terms  proportional  to 
gives  the  following  differential  equation  for  /„  (y  ,5) : 


<9/ 


(rl + Mo<^nn)L  =  0  • 


The  solution  for  /  is 

*'  n 


(43) 


(44) 


where 


K=^|Yl  +  ^l■o<^s+jf^,ovy„  . 


(45) 


Since  the  magnetic  field  must  go  to  zero  at  y-oo,  we  require  that  E„=0  and  thus 
have 


Inserting  equation  (46)  into  equation  (41)  gives 


n= 


(46) 


(47) 


Using  the  same  methodology  for  equation  (29b),  we  derive  the  following  expression  for 
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(48) 


^  n=+oo^  * 

By=  I 

n=-oo 

where  ^(5)  is  related  to  ^(5)  through  the  divergence  relationship  of  equation  (30): 

dB  n--K>o 

^+-^=  I  (jr.F.(s)-KP,y‘”'e>'-’  =0.  (49) 

The  only  way  that  equation  (49)  can  be  satisfied  is  if  the  “(  )”  terms  are  equal  to  zero. 

Setting  these  terms  equal  to  zero  gives 

4  =  (50) 

From  the  foregoing  equations,  we  see  that  the  F„ ’s  determine  the  x-  and  y-components  of 
the  electromagnetic  field.  The  F„ ’s  are  determined  fi'om  the  boundary  condition  of 
equation  (36)  applied  in  Laplace  transform  space.  Recalling  that 

15(0  =  1,  (51) 

where  L  is  the  Laplace  transform  operator,  and  evaluating  equation  (47)  at  y  =  0  gives 

n~oo  Ort'TTY 

s: ■  (52) 

n=l  ^  n=-<*> 

Using  =cosy„A:+  7sin7„A:  gives 

^  =  ^  [(/;  +  F.„ )cosy„x  +  j{P„  -  F_„ )sinr„x] .  (53) 
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Inserting  equation  (53)  into  equation  (52)  and  noting  that  there  is  no  Qq  term  in 
equation  (52)  gives 


+  (54a) 

F„-F_„=0,  (54b) 


and 


^0=0. 


(54c) 


From  equations  (54a  -  c),  we  have 


(55) 


The  time  domain  A:-component  of  the  magnetic  field  is  obtained  by  taking  the  inverse 
Laplace  transform  of  equation  (47): 

\  F„it-r)Q„{r)dte^"\  (56) 

n=-oo  ^ 

where  F^(t)  is  the  inverse  Laplace  transform  of  F„(s)  and 

G„(T)  =  =  F-\(s) .  (57) 

Qfi(x)  is  computed  in  two  steps.  First,  we  write 
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where 


+>r„- 


(59) 


Then 


(60) 


From  the  properties  of  the  inverse  Laplace  transform,  we  have 


(61) 


where 


(62) 


Using  the  formula 


2^  TV  t 


then  yields 


a.(T)=e- 


e-'-'>Y(y,T) 


:y  >0. 


(64) 


Equation  (56)  is  simplified  using  the  result 
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F„(t)  =  F_„(t)  = 


2 


S(t). 


(65) 


Substituting  equations  (64)  and  (65)  into  equation  (56)  yields 


n=+<» 

s 


«=-oo 


Q„|  5{t-T:)e ''’'‘'\jr{y,x)e^’’‘dT . 

o 


(66) 


Using  the  formulas 


- 


(67a) 


r-n  = 


y« 


■Prn  ='■«» 


(67b) 


and 


r-„  =-r«  = 


—2nK 

"T”’ 


(67c) 


reduces  equation  (66)  to  the  simpler  forms 


r;T 


S(t-Tyiriy,r)e  cos(y„x-y„vT)dr 

n=i  o 


(68) 


and 


7l* 


W(yd)e  cos(7„x - 7„vr) 


n=l 


(69) 
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We  now  compute  By(x,y,t)  by  taking  the  inverse  Laplace  transform  of  equation  (48). 
Using  equations  (50),  (54a  -  c),  and  (58),  we  initially  have 


„“L  4B^yls  +  r„  [  2 


P'o^n  1 


followed  by 


S(t-T)^  y,Q,e^-T.(z)dT . 

^  V  O'  „  n=-«> 


where 


r„(T)  =  L-'  1 

Jj  +  r. 


Using  equation  (61),  we  get 


-Mr 

r„(T)  = 


-.r  1 


4t  _ 


e  '■"*$(y,T) , 


where 
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I 

<I)(y,T)=-| — e 

Jtct 


Substituting  equations  (47),  (73),  and  (74)  into  equation  (71)  and  then  integrating  over  t 
gives 


_  ;  III  1 

-rJ-T=*  *' 

2  \  a 


Since  y„  is  an  odd  function  of  n,  while  is  and  even  function  of  n,  we  get 


Combining  positive  and  negative  values  of  n,  and  using  equation  (76)  gives 


+  y.  A,^"""'"”’  =  2#A  sin(y,A:  -  r.vO  • 


,+j7-.(x-vt)  _ 


Inserting  equation  (77)  into  equation  (75)  yields  the  following  expression  for  B  : 


By  =  ^  ^  Sin(y„x  -  y„vt) 


ujjy  -M. 


Lastly,  the  current  density  in  the  moving  conductor  is  determined  from  the  formula 


j  =  VxH  =—VxB. 
1^0 
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4.  Characterization  of  External  Magnetic  Field 


The  purpose  of  this  section  is  to  calculate  the  environmental  magnetic  field  in  the 
region  below  the  source  conductors  (see  Figure  1),  that  is,  for  field  points  satisfying  the 
condition  yo<-w.  The  magnetic  field  at  the  point  XQ,yQZo  is  determined  in  a 

straightforward  way  using  the  Biot-Savart  formula  of  equation  (1)  with  x,y,z  replaced 
by  jcq,  Jo  Zq  .  The  volumetric  integration  in  equation  (1)  is  over  all  current  densities. 

In  our  problem,  there  are  two  current  sources  that  contribute  to  the  environmental 
magnetic  field:  the  first  is  the  source  firom  the  input  current  source  elements  and  the 
second  is  from  the  spatial  distribution  of  current  density  within  the  moving  conductor. 

Using  the  formula 


VxB  =  it 


dx  dy 


(80) 


the  z-component  of  current  density  in  the  moving  conductor  is  given  by 


Jc=- 


dB^  dB, 


fi„[dx  By 


(81) 


The  z-component  is  the  only  component  of  current  density  in  our  system  and  is  a 
consequence  of  the  geometry.  Using  the  results  of  section  3,  which  provide  B^  and  By  as 

a  function  of  x,  y,  and  r,  we  can  calculate  the  spatial  distribution  of  current  density  for  the 
moving  conductor.  When  the  resulting  expression  for  j^is  inserted  in  equation  (1),  the 
environmental  magnetic  field  from  the  moving  conductor  is  determined  as  a  function  of 
Xq-,  yo»  ^0’  3nd  r  in  the  region  yg  <->»'- 


27 


In  contrast  to  the  exact  computation  of  the  environmental  field  from  the  spatial 
distribution  of  current  density  within  the  moving  conductor,  the  contribution  from  the 
source  current  is  relatively  easy.  The  y-dimension  of  the  source  current  is  assumed  to  be 
very  thin  compared  with  |yo  ”  H  •  moving  conductor,  the  vertical  integrated 

current  (i.e.,  image)  is 


J*=]jc{y)dy=^ 

0  “o 


(82) 


Placing  /*  at  some  equivalent  height  above  the  ground  plane  y  =  0  is  analogous  to  placing 
an  image  current  a  distance  w^=w  +  d  above  a  finite  conducting  ground  plane. 

From  the  infinite  conductivity  case,  we  know  that  the  current  in  the  conductor  will  act 
to  compensate  the  magnetic  field  generated  by  the  source  current.  For  this  a  =  oo  case, 

O 

the  integrated  vertical  current  is  a  continuous  planar  distribution,  180  out  of  phase  with 
the  source  current.  This  image  is  located  a  distance  w  above  the  y  =  0  plane.  The 
resultant  environmental  field  will  be  the  smallest  in  this  case  because  it  is  composed  of 
two  parts  that  nearly  cancel  each  other  when  |yo|  <  -w . 

When  <7  oo ,  the  image  current  cannot  be  placed  any  closer  to  y  =  0  than  w.  This  is 
true  irrespective  of  the  velocity.  In  the  other  extreme,  we  can  place  the  image  current  so 
far  above  the  y  =  0  plane  (this  is  only  a  mathematical  abstraction)  that  the  environmental 
magnetic  field  will  be  due  to  only  the  source  current  elements. 

hi  summary,  the  upper  bound  on  the  environmental  magnetic  field  will  be  determined 
by  completely  neglecting  the  current  in  the  moving  conductor.  The  lower  bound  will  be 
computed  by  locating  the  image  current  a  distance  w  above  the  y  =  0  plane.  We  choose 
to  render  the  results  in  this  format,  not  because  we  cannot  calculate  the  distant  magnetic 
fields  exactly  [^can  be  calculated  exactly  from  equation  (81)  for  use  in  the  Biot-Savart 
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law],  but  because  this  viewpoint  provides  enhanced  insight  into  the  salient  physical 
processes. 


Equation  (82)  can  be  formally  integrated  to  yield 

/*  =  —\G{x,t)-B^{x,y  =  oo,t)  +  B^(x,y  =  0,o], 


(83) 


where 


dBy(x,y,t) 

dx 


dy . 


(84) 


Taking  the  partial  of  equation  (78)  with  respect  to  x  and  then  integrating  over  y  gives 


Gix,t)=—  cos(y„x-r„vr) .  (85) 

^  n-l 


The  exponential  terms  in  equation  (85)  are  a  result  of  diffusion  in  the  x-direction.  The 
modal  diffusion  time  constants  are  given  by 


4n^7C^ 


(86) 


Using  iig  =  Any.  10'^ Him,  o  =2.2yl(f  mhosim,  and  a  nominal  dimension  of  L  =  0.5  m 
gives 


tji  =  0.18  s 


(87) 
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for  the  first  mode.  Even  though  higher-order  modes  will  have  shorter  diffusion  times, 
they  will  generally  be  much  longer  than  the  pulse  lengths  of  interest — ^typically  in  the 
milUsecond  range.  We  return  to  this  time-constant  issue  after  our  general  analysis  is 
completed. 

Returning  to  equation  (83),  we  note  that  since  B^(x,y  =  °o,t)=  0,  we  can  use 
equations  (34)  and  (35)  to  write 


Gjxj) 


-5(0S^«cosy„jc. 


n=l 


(88) 


Because  all  the  magnetic  field  components  determined  up  to  this  point  were  based  on 
the  impulse  function,  so  is  /* .  Applying  the  integral  relationship  of  equation  (37)  as 
applied  to  J* ,  the  step  function  response  for  the  total  current  per  unit  length  is 

t  t 

Cep  =  aJ/*(/  -  T)dT  =  A  J  J\x)dT .  (89) 

0  0 


When  equation  (88)  is  inserted  into  equation  (89),  we  get  two  contributions,  which 
are  distinctly  associated  with  the  image  contribution  described  in  section  2  and  another 
part  that  is  identified  with  diffusion  and  motion.  We  write 


7  =J  -J 

''step  "'dm  •'/  » 


where  the  image  contribution  is  given  by 


/  n=oo  n=eo 

/;  =  aJ5(t)JtX^„  cosy„x  =  A^X^n^osy^x, 

0  «=1 


(90) 


(91) 


and  the  diffusion-motion  contribution  is  given  by 
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(92) 


JlM=^]G{x,x)dT. 

r^o  0 

Equation  (92)  is  evaluated  using  equation  (85).  Before  presenting  the  general  result,  it 
is  instructive  to  evaluate  for  the  special  case  of  a  stationary  conductor,  v  =  0.  These 
results  confirm  the  assertions  rendered  in  section  2.  When  v  =  0,  we  have 


G(x,t)=  — 
a 


n=+<»  — YsL 

‘^ylQ.„cosy„x, 


(93) 


and 


A  r  S'  — 

C  =  — I  G{x,x)dT  =  A,  (1  -  e  '^ )  cos7„x . 

0  "=i 


(94) 


In  carrying  out  the  integration  in  equation  (94),  we  have  used  the  formula 


(95) 


Equations  (91)  and  (94)  are  rich  in  physical  content  and  can  be  traced  back  to  the 
discussion  of  section  2.  If  we  let  tr  =  <» ,  then  ~  X  )”  t®>nns  of 

equation  (94)  vanish,  and  we  are  left  with  only  the  image  contribution, 


(96) 


By  comparing  equation  (94)  with  equation  (91)  on  a  term-by-term  basis  and  noting 
that 
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_ ^ 

(97) 

we  readily  show  that,  on  the  average,  over  the  periodic  spatial  interval  {0,  L},  the 
magnitude  of  is  less  than  the  magnitude  7; .  The  proof  is  straightforward  using  the 
orthogonality  of  the  cosine  terms.  We  have 

,  (98) 

2  ) 


and 


For  all  real  physical  systems,  J]  must  be  bounded  and,  hence. 


(99) 


(100) 


Every  term  in  the  series  of  equations  (98)  and  (99)  is  positive.  Since  each  term  in 
equation  (99)  is  greater  than  the  corresponding  term  in  equation  (98),  the  proof  that 

kf  I  >  is  complete.  In  the  regime  where  t«t.,  we  have  the  approximation 

I  ^  \ave  ' 


^tlt^  «l. 


(101) 


In  this  case,  image  theory  provides  a  very  good  description  of  current  flow. 
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Lastly,  we  notice  that  equations  (91)  and  (94)  can  be  inserted  in  equation  (90)  to  give 
the  result  for  a  stationary  conductor  based  on  a  first-order  perturbation  to  the  perfect 
conductor  solution  aXy  =  0 


'^step  ^ 

/t=5l 


cosy„A: . 


(102) 


Equation  (102)  shows  that  the  current  in  the  conductor  will  eventually  vanish  for  a  step 
function  excitation  and  is  due  to  conductive  energy  loss. 

We  now  consider  the  general  case  involving  a  moving  conductor,  which  requires  use 
of  equation  (85)  in  equation  (92).  This  calculation  is  facilitated  by  first  writing 

cos(y„x  -  Y„vt)  =  cos(y„x)  cos(7„vr)  +  sin(7„x)  sin(y„v0 .  (103) 

Inserting  equation  (103)  into  equation  (85)  gives 

G(x,t)  =  Gi(x,t)  +  G2(x,t),  (104) 


where 


■t  n=4<o  . 

G,U.0  =  -S« 

^  »ssl 


J2. 


r„a„cos(r„x)cos(y„v0 


(105) 


and 


11=4^  — 5j5f_ 

Sin(y„x) sin(y„v0 . 


(106) 


Substituting  equations  (104)  -  (106)  into  equation  (92)  and  integrating  yields 
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(107) 


JIm  =  COS(r„J^)¥„,(0  +  -^Yr^n^n  Sin(Y„X)W^it)  , 

«==!  «=I 


4 


where 


1 


''  rnt 

f  y^  V 

HcO 

-^cos(7„vr)  -  7„vsin(7„v0 

Jj 

(108) 


and 


1 


-^f  -v2 


y„v-e 


^sm(j„vt)  +  y„vcos(r„vt) 


.  (109) 


The  foregoing  expressions  show  that  the  conductor  current  contains  harmonics  with 
angular  frequencies: 


«„=r„v=- 


InTtv 


(110) 


Using  values  of  L  =  0.5  m,  and  v  =  200  m/s  gives 

0)i  =2,500  radls 


(111) 


for  the  first  harmonic. 

The  damping  time  constants  given  by  equation  (86)  show  that  the  higher  harmonics 
dampen  out  much  more  rapidly  than  the  lower  numbered  ones.  This  damping  produces 
frequency  broadening,  which  must  be  taken  into  account  when  anal5^ing  the  response  of 
electronic  equipment  that  is  located  near  the  system  of  conductors  shown  in  Figure  1. 
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In  order  to  better  understand  the  essential  features  of  environmental  magnetic  fields 
generated  by  the  moving  conductor  current  system  of  Figure  1,  it  is  desirable  to  compare 
the  magnitude  of  the  modal  terms  in  equation  (107)  with  the  corresponding  modal  terms 
of  J]  as  given  by  equation  (91).  We  shall  use  the  same  procedure  for  making  this 
comparison  as  we  did  in  the  development  of  equations  (98)  and  (99). 

Using  the  orthogonality  of  both  cosine  and  sine  function  in  computing  froiti 

equation  (107),  we  obtain 


2  n=«> 
^  n=l 


\i;2 


(112) 


where 


H 


2 

n 


ns 

^dn 


(113) 


With  the  help  of  basic  trigonometry,  algebra,  and  equations  (108)  and  (109),  we 
readily  prove  that 


H  =■ 


n 


(114) 


where 


IL  _L 

h^(t)  =  \  +  e  -  2e  **  cosa)„t . 


(115) 
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Figures  5  and  6  show  A„(Oand  H^(0,  respectively,  as  a  function  of  time  with  mode 
number,  n,  as  a  parameter  for  =  47txl0'^ Him,  <7  =2.2xl(f  mhos/m,  L  =  0.5  m,  and 
V  =  200  m/s. 


Figure  5.  h„  (t)  as  a  Function  of  Time  With  Mode  Index,  n,  as  Parameter. 


Figure  6.  H^(r)  as  a  Function  of  Time  With  Mode  Index,  n,  as  Parameter. 
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By  comparing  equation  (112)  with  equation  (99)  on  a  term-by-term  basis,  we  note  that 
their  ratios  are  H^.  Equation  (115)  shows  that  the  maximum  value  of  \it)  is  four. 

Thus,  for  this  exploratory  study,  with  the  mode  index,  n,  less  than  10,  H^can  be 
approximated  by  its  upper  bound  as 


H 


2 

n 


4n^ 

n  +(t^j6)j) 


(116) 


We  also  have 


=  0.18  (2,500)  =  450. 


(117) 


The  higher  order  terms  of  1/^^  I  become  comparable  to  those  of 


when  H, 


approaches  unity.  Equations  (116)  and  (117)  show  that  this  will  occur  when  the  mode 
index  becomes  greater  than 


n  = 


1/2 


=  260. 


(118) 


Considering  the  first  few  modes,  |7^j^|^^will  be  significantly  less  than  The 

motion  of  the  conductor  does  not  appreciatively  affect  the  environmental  field. 


For  the  purposes  of  calculating  and  bounding  the  environmental  magnetic  field  in  the 
region  <  —w  we  need  to  consider  the  real  current  density  source  located  in  the  plane 
y  =  -w  and  its  image  current  density  source  located  in  the  plane  y  =  +w  (see  Figures  1 
and  2).  As  previously  mentioned,  the  maximum  environmental  magnetic  field  will  occur 
when  the  image  contribution  is  absent  and  the  smallest  will  occur  when  the  image  is 
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included.  These  two  limits  bound  the  problem,  irrespective  of  how  well  we  calculate  the 
details  of  the  current  in  the  moving  conductor. 

The  starting  point  for  the  calculation  of  the  magnetic  field  at  position 

r„=^x„  +  ]y„  +  kz„  (119) 


is  the  Biot-Savart  law,  which  is  given  by  equation  (1).  In  order  to  be  consistent  with  the 
notation  of  section  3,  we  let  x,  y,  z  be  the  coordinates  of  the  volumetric  distribution  of 
current  density,  the  reader  may  recall  that,  in  section  2,  x,y,z  were  used  for  the 
volumetric  distribution  of  current  density.  When  we  consider  only  the  real  source  and  its 
image,  the  magnetic  field  at  is  given  by 

=  +  (120) 


where  Bj^  and  Bj  are  the  contributions  from  the  real  and  image  sources,  respectively. 
They  are  given  by 


■+00  +h 


-oo  -A 


(121) 


and 


4;r 


(122) 


-A 


r/ 


where  the  integration  over  x  extends  from  - to  +  <» ,  the  integration  over  z  extends 
from  -h  to  +h: 
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rR=T{x^-x)  +  ~j{,y,  +  w)  +  k (z„  - z) , 

(123a) 

*  X =  j(Xo-x)-i(yo  +  w) , 

(123b) 

=  ((Xo  -xf+(yo  +  +  (4  -  z?  , 

(123c) 

n=T{x^-x)  +  jiy„-w)  +  kiZo-z), 

(123d) 

X 

II 

N 

1 

T 

1 

(123e) 

1^/  f  =  ((^«  -xf+iy„-  wf  +  (z„  -  zf  . 

(123f) 

In  equations  (123a  -  f),  we  place  no  restriction  on  the  location  of  the  field  point,  although 
we  apply  the  results  to  those  locations  where  y,,  <  -w . 

The  negative  sign  in  equation  (122)  and  the  appearance  of  in  the  integral  accounts 
for  the  fact  that  the  image  current  is  the  negative  of  the  real  current.  Equations  (123a  -  f) 
accomit  for  the  fact  (see  Figure  2)  that  the  source  current  density  is  located  at  y  =  —w , 
while  the  image  current  density  is  located  at  y  =  +w.  Jj^  is  determined  from 
equations  (21)  and  (24),  with  given  by  the  step  function  of  equation  (38).  We  have 

'2,nltx 

•^«W  =  ZA.cos— -  (124) 

n=l  ^ 


and 


(125) 


2A, ,  .  .  nnh 

A,  =  — =^(1  -  cosn;r)sin— — . 

^  nit  L 

The  calculation  of  the  analytical  forms  for  and  Bj  is  straightforward.  For  brevity, 

we  calculate  and  then  simply  determine  Bj  by  substituting  y^-w  in  place  of 
+  wand  taking  the  negative  of  the  result.  Equation  (121)  is  first  integrated  over  z  to 
yield 


=  -I  {yo  +  ’ 


(126) 


where 


+h 

w^U 


dz 


L[b^  +  {z„-zf) 


3/2  » 


(127) 


and 


b'^  +{y„  +  wf . 

The  integration  in  equation  (126)  gives 


(128) 


Zp-^b 


z„-h 


((z,  +  hf  +  "  ((z,  -hf+b^  y " 


(129) 


The  foregoing  expression  simplifies  when  we  take  our  observation  point  to  be  located 
at  the  midpoint,  z^  =  0.  There  is  essentially  no  loss  of  generality  in  this  case,  and  W 

becomes 
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h 


(130) 


(h^+bT 


Since  one  of  the  principle  aims  of  this  study  is  to  examine  the  relationship  between 
the  distance  of  the  field  point  from  the  conductors  as  compared  to  the  spatial  periodicity 
of  the  current  in  the  moving  conductor,  we  further  simplify  W  by  setting  =  0 . 

Equation  (130)  then  reduces  to 


W  = 


2h _ 1 _ 

+  (yo + [h^ +x^  +  (yo  + 


Equation  (131)  is  used  in  subsequent  calculations. 


(131) 


The  procedure  from  this  point  is  to  substitute  equations  (124),  (126),  and  (131)  into 
equation  (121)  and  then  integrate  over  x  between  the  limits  of  —  <»  to  +<» .  Formally 
integrating  over  x  initially  yields 


B,= 


4k 


+0O 

'  iyo  +  >^)  J  + J  J  JR^wdx 


(132) 


Since  Jg  and  W  are  symmetric  functions  of  x,  and  the  “x”  in  the  second  integral  is  an 
antisymmetric  function  of  x,  the  second  integral  vanishes.  Inserting  equations  (124)  and 
(131)  into  equation  (132)  gives 


R  -  IL 


i(yo + w)S4®«(>o  + 


n-l 


(133) 


where 
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t-oa 

®,  =  J 


Innx 


2h 


cos- 


L  x^+  +  wf  +  x'^  +  (y^  +  wf 


dx. 


(134) 


Simplification  results  when  we  cast  the  foregoing  equations  in  dimensionless  form  by 
introducing  the  variables 


P 


(135) 


77  = 


yo+^ 

L 


(136) 


and  selecting  a  specific  value  for  h.  For  convenience  a  value  of  h  =  Lis  used. 

Using  equations  (125),  (135),  and  (136)  we  derive  the  following  expression  for  , 
recalling  from  equation  (133)  that  it  is  the  x-component  (hence  we  drop  the  vector 
notation): 


(137) 


where 


a  =  —(1  -  cosn;r)sin-^^^ 

"  HTC  L 


(138) 


and 


coslnjtp 

{jf+p^)(l  +  T}^  + 


dp. 


(139) 
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Figure  7  shows  as  a  function  of  positive  77 ,  with  mode  index,  n,  as  parameter.  As 
expected, decreases  with  increasing  distance  from  the  moving  conductor  and  with 
increasing  mode  index.  We  also  have 

=  (140) 

which  is  relevant  in  our  calculations  since  is  negative. 

5 
4 
3 

Pn 

2 
1 
0 

0  0.2  0.4  0.6  0.8 

F%ure  7.  P„(tj)  as  a  Function  of  Positive  77  With  Mode  Index  as  Parameter. 

The  total  magnetic  field  is  computed  from  the  formalism  of  equations  (120)  to  (123). 
In  terms  of  the  field  point  used  in  the  foregoing  calculations  =  0 ,  =  0 ,  we  have 

(141) 

the  second  term  being  attributed  to  the  image. 
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5.  Numerical  Results 


The  purpose  of  this  section  is  to  obtain  numerical  estimates  of  the  magnetic  flux 
density  as  a  function  of  vertical  distance  from  the  source  conductor  system  for  a  particular 
set  of  parameters.  The  following  values  are  used: 

/,  =  175  mm. 


I2  =  75  mm, 

L  =  2(/,+/2)=  0.5  m, 

w  =  2.5  mm, 

h  =  L, 


and 


4  =  600,000  A. 

There  is  a  wide  range  of  magnetic  field  variations  in  this  system.  The  highest  fields 
occur  very  close  to  the  source  elements  themselves.  In  this  region,  which  is  not 
anticipated  to  be  close  to  any  external  electronics  systems  in  practical  applications,  the 
magnetic  field  intensity  is  of  the  order  of 

H  ^  —  =  A,  =  3.4  kA/mm  (142) 

/, 


and  the  magnetic  flux  density  is 
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B  =  =  4.2  T. 


(143) 


As  we  move  away  from  the  individual  conductors,  the  influence  of  the  periodic  array 
of  the  source  ciuxent  is  seen  and  the  theory  provided  in  the  previous  section  must  be  used. 
The  largest  fields  in  this  case  exist  when  the  image  effect  is  negligible.  This  occurs  when 
the  image  distance  becomes  large.  As  previously  discussed,  the  case  in  which  the  image 
currents  can  be  neglected  is  a  mathematical  abstraction  and  provides  a  limiting  case  of 
maximum  environmental  magnetic  field.  This  field  is  given  by 

-^AjYanPn(J]) .  (144) 

Using  the  numbers  provided  at  the  beginning  of  this  section  and  noting  from 
equation  (138)  that 


=  0  :  n  =  even  number. 


(145) 


we  obtain  a;  =  1.15  and  a;  =  -0.064. 

Subsequent  values  of  are  negligible.  Moreover,  when  aj  and  as  are  combined  with  the 
results  of  Figure  7,  we  see  that  the  third  harmonic  is  just  a  few  percent  of  the  first 
harmonic.  Hence,  with  very  good  approximation,  we  can  use  the  result 


Alt 


(146) 


where  the  dimensionless  distance  is 


rj  = 


(147) 
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Figure  8.  Magnetic  Field  as  a  Function  of  Dimensionless  Distance. 

In  the  foregoing  equation,  we  have  adjusted  for  the  fact  that  is  negative.  Figure  8  is  a 
plot  of  (without  the  image  current)  as  a  function  of  ij . 

When  the  image  is  taken  into  account,  we  must  use  equation  (141),  which  now 
becomes 


(148) 

A  plot  of  equation  (148)  is  also  shown  in  Figure  8.  As  observed,  the  magnetic  field  with 
the  contribution  from  the  image  current  reduces  the  environmental  field  by  more  than  an 
order  of  magnitude.  However,  in  both  cases,  the  field  has  decayed  roughly  3  orders  of 
magnitude  at  one  conductor  spacing  (L)  from  the  source. 
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6.  Summary  and  Conclusions 


We  have  developed  a  model  for  determining  the  space-and-time  behavior  of  the 
magnetic  field  external  to  a  system  of  moving  pulsed  conductors.  The  model  includes  a 
moving  finite  conductivity  conductor  in  relative  close  proximity  to  a  stationary  conductor. 
The  analysis  indicates  that  the  magnetic  field  has  a  space-and-time  periodicity  that  can  be 
related  to  the  parameters  of  the  conductors.  Furthermore,  we  present  a  numerical 
example  that  predicts  the  magnitude  and  frequency  of  the  magnetic  field  as  a  function  of 
distance  from  the  stationary  conductor.  In  this  analysis  we  show  that  the  induced  currents 
are  nearly  equal  in  magnitude  to  those  in  the  source.  The  eddy  currents  in  the  moving 
conductor  delay  the  penetration  of  the  magnetic  flux.  The  environmental  magnetic  field 
decays  rapidly  to  negligible  levels  at  distances  of  interest  to  the  EMC  problem. 

One  of  the  more  interesting  and  potentially  important  devices  currently  in 
development  is  the  air-core  pulsed  alternator.  These  devices  have  been  used  to  deliver 
large  currents  for  electromagnetic  launch  systems.  The  main  difference  between  a 
conventional  generator  and  a  pulsed  alternator  is  die  latter’s  extensive  use  of 
compensation.  Compensation  is  a  method  of  reducing  the  inductance  associated  with  an 
alternator’s  armature  winding.  This  reduction  is  achieved  by  inducing  currents  in  a 
moving  conductor  that  is  in  close  proximity  to  the  armature  winding  of  the  alternator. 

hi  order  to  estimate  the  environmental  magnetic  field  produced  by  a  moving  pulsed 
conductor,  we  have  assumed  a  step  function  current  for  the  current  source  elements,  hi 
an  actual  alternator  system,  these  currents  would  be  produced  as  a  result  of  voltage 
generation  due  to  a  spatially  varying  direct  current  excitation  field  established  in  the  air 
gap  between  the  moving  conductor  and  the  armature  windings.  Of  course,  this 
additional  source  would  contribute  to  the  environmental  magnetic  field.  This  excitation 
source  has  not  been  considered  in  this  analysis  but  can  readily  be  computed  using  the 
techniques  developed  in  this  study. 
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